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Abstract 

Larval zebrafish exhibit a variety of complex undulatory swimming patterns. This repertoire is controlled by 
the 300 neurons projecting fi-om brain into spinal cord. Understanding how descending control signals shape 
the output of spinal circuits, however, is nontrivial. We have therefore developed a segmental oscillator model 
(using NEURON) to investigate this system. We found that adjusting the strength of NMDA and glycinergic 
synapses enabled the generation of oscillation (tail-beat) frequencies over the range exhibited in different larval 
swim patterns. In addition, we developed a kinematic model to visualize the more complex axial bending patterns 
used during prey capture. 
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1. Introduction 

The spinal cords of vertebrate animals contain segmental oscillators, or central pattern generators 
(CPGs), that can produce rhythmic movements. The operations of these spinal CPGs are best under- 
stood in l ower vertebrates , such as lamprey an d Xenopus, where they are used for undulatory swimming 
|)R obert j . Il998l: l&iichananL Il999l: Ic^illneil Efloil) . Control signals descending from brainstem to spinal cord 
have also been studied extensively in these and other lower vertebrates, such as goldfish and zebrafish. 
Their lack of a corticospinal tract avoids a degree of complexity that is present in mammals. The func- 
tioning of descending control systems in higher vertebrates has been difficult to understand. Studies of 
lower vertebrates should reveal conserved principles by which these systems operate. The larval zebrafish 
takes vertebrate simplicity to an extreme: the decreased numbers of neurons allows exact identification of 
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many cell types in both brainstem and spinal cord l|Kimmel et al l llflSlllflS.'llLiu and WesterfieldL llQsi 
iBernhardt et allll990() : this in turn provides major experimental and modeling advantages. 

One important aspect of lower vertebrate locomotion is the regulation of swimming speed, which is often 
correlated with the frequency of alternating left and right contractions of the axial muscles of the trunk or 
tail, termed tail-beat frequency (TBF). In a steady-swimming fish TBF is generally equal to the oscillation 
frequency of the spinal CPGs, so by understanding the modulation of CPG frequency in spinal cord we 
can understand a major element of the control of swim speed. In goldfish, lamprey, and other fishes, 
swim frequency can be modulated by stimulation of the midbrain locomotor region or bath application 
of NMDA; It is also known that serotonin, acetylcholine, dopamine, and other neurotransmitters also 
influence the CPG's oscillation frequency. However, the identities and locations (in brainstem or spinal 
cord) of the cells involved in swim frequency control are unknown. Larval zebrafish may provide insights 
into this problem because they exhibi t distinct swim patterns tha t span a broad range of tail beat 
frequencies, ranging from 25 to 75 Hz ijBudick and O'MallevL fioOCf) . A critical unknown is whether or 
not the distinct larval swim patte rns (slow, bu rst, and cap ture) are c ontrolled by distinct control systems 
and/or specific neurotransmitters l|Borla et al.t . ,2Q02.; O'Mallev et al.i^.20Q3) . Modeling the different types 
of swimming behavior using a combined neural and kinematic model can shed light onto the presence 
and relevance of different motor control elements. 

The functionin g of spinal networks that underlie locomo t ion in fishes and tadp oles has been extensively 
modeled (see e.g. iRoberts and TunstallL Il99(it IPalel ll995L l2003l iGrillneil 1200,'^ . We created a zebrafish 
neural model, based on previous Xenopus spinal network models, of the CPGs i n the larval z ebrafish spinal 
cord, incorporating known properties of the oscil l ators underlying swimming IjTunstall et al.. 2002). The 
spinal interne uron types in zebrafish ( Hale et all l200ll) are likely homologous to those in both Xenopus 
and lamprey i Fetchol Il992(l . We explored the control of tail-beat frequency, and found that by altering 
the strengths of AMPA, NMDA and glycinergic-like synapses (all known to be present in lower vertebrate 
spinal cords), we were able to generate TBFs that spanned the range of speeds observed during burst 
and slow swimming behaviors. We also created a simple mechanical or "kinematic" model (that can be 
driven by the neural model) to visualize how the spinal neural activity might be transformed into larval 
behaviors. Our ultimate goal in creating this neurokinematic model is to introduce a tool for testing 
theories of descending motor control in the larval zebrafish. 



2. Methods 



Figu re 1 shows the structure of our neural model, which is a simplification of the larval spinal cir- 
cuitry 1 Hale et all 12001) incorporating the minimal elements required to generate rhythmic, propagating 
alternating activity. Each segmental oscillator consists of two neurons, with NMDA and AMPA-like au- 
tosynapses, which are connected by glycinergic synapses. An oscillator can be triggered with a single 
excitatory pulse to one cell, followed by a single pulse to the other cell several milliseconds later (7 ms 
in our simulations). Individual oscillators are connected into a 25-segment chain (corresponding to the 
approximately 25 segments in the zebrafish spinal cord), with nearest- neighbor descending excitatory 
synaptic connections; the entire chain can be started merely by triggering the head segment. 

We use cells with the standard Hodgkin-Huxley channels provided by the NEURON soitwaie package: 
persistent potassium (gx = 0.036 S/cm^, Ek = — 77mV), transient sodium (gNa = 0.12 S/cm^, i?Na = 
50 mV), and leak {gj^ = 0.0003 S/cm^, i?L — — 54.3mV) channels. Our synap ses are modeled using the 
difference of two exponentials, with rise time ti(= 1 ms) and fall time T2. See lDavan and Abbott! l|200lL 
p. 182) for details. The three types of synapses we used were AMPA (t2 = 6 ms, E = OmV), NMDA 
(t2 = 80 ms, E = 60 mV), and glycine (t2 =2ma, E = -80 mV). 
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Fig. 1. Schematic diagram of the neural oscillator model. The squares represent neurons, and each line indicates one of three 
types of synaptic connection: glycinergic (fast inhibitory), AMPA-like (fast excitatory), or NMDA-like (slow excitatory). 
Only the first two segments of a 25-segment chain are shown. 

Our kinematic model is meant to translate the neural signals from the spinal circuitry into observable 
kinematic behaviors. To do so, we use a relatively simple transformation of putative neural output into 
effects on the radius of curvature of a line segment representing the trunk of the larva. This model does 
not take into accoimt the full physics (e.g., elasticity and hydrodynamics) of the situation; nevertheless, 
it does approximate the observed axial kinematics recorded experimentally with a high-speed camera. 

Suppose that each segment x is receiving a neural signal Fs{x,t) at time t from the spinal circuitry: 
Fs > for a signal to the right side of the segment, Fs < for a signal to the left. We suppose that this 
signal is integrated through exponential synapses, so that the signal passing to the muscles is 



r 

Fm{x,t)= J F,{x,t')[e^''- 
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(1) 



where n and T2 are the growth and decay time constants of the synapse (we use n 
in our calculations.) Assuming that the muscle contracts linearly as a function of F, 
radius of curvature of segment x is 

Rix,t) = Wix)/Fmix,t) 



6 ms and T2 = 8 ms 
m, it follows that the 



(2) 



W{x) is a function describing the stiffness, or resistance to bending, of segment x. It takes into account 
the width of the body (the tail is more flexible, and bends more than the rostral trunk). 

We can feed the output of our neural model directly into this system, but it is also useful to introduce 
a (more) artiflcial signal. We build this signal out of three components: an oscillatory signal Fosc{x,t), 
which is a series of delta functions propagating caudally; a bending signal Fjend, which is a tonic signal 
applied to one side of all segments at once; and a rostral stiffening signal, which reduces the signal to 
the first Xinh rostral segments by a factor finh- The first two pieces are responsible for swimming and 
turning, respectively. The stiffening signal is seen in larvae during prey capture, where the fish keeps its 
head relatively still while adjusting its orientation with its tail. We can write our artificial neural signal 
as 



Fs{x,t) = {Fosc{x,t) + fbend) X 



J finhi 



X !> Xinh 
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where fbend is a constant and 



3 



N 
I 



3 
(X 
0) 



70 
65 
60 
55 
50 
45 
40 
35 















9nmd 


, = 0. 


B3 nS 
















9nmd 


, = 0. 


40 nS 




















WMMMWI 














9nmd 


, = 0. 


25 nS 
















9nmd 


, = 0. 


16nS 


























1 1 


1 1 


I 1 




1 1 


~I 

1 1 


1 


1 1 


— •■- 
1 1 





100 200 300 400 500 600 700 800 900 1000 

time (ms) 

Fig. 2. Instantaneous TBF in a single segmental oscillator, for several values of gNMDA- Each point represents the iTBF for 
the left neuron; the corresponding points from the right side gives near- identical results. Oscillators approach a steady-state 
frequency only after an initial "wind-up" period. 



1, X = lnvt (mod A) 
Fosc{x, t) = fosc "( -1, X = 2Tivt -I- A/2 (mod A) 
0, otherwise 



(4) 



The parameter v is the tail-beat frequency, and A is the length (in number of segments) of the resulting 
wave which propagates down the fish. 



3. Results and Discussion 



We began by studying a single segment of the neural model shown in Fig. 1. A transient trigger 
pulse to this model initiates sustained alternating activi ty, mimicking Xenopus spin al cord, where a 
transient stimulus produces a sustained bout of swimming feobert s and TimstallLllQQ O^. The first larval- 
zebrafish specific task was to generate the range of tail-beat frequencies (TBFs) used in different swimming 
behavior s. The frequency o f oscillation was calculated on a half-cycle basis (termed "instantaneous" 
TBF in iBorla et alJ.l2002(l . Figure 2 shows how instantaneous TBF varies with time, with each trace 
representing the frequency profile for a specific value of the NMDA synaptic conductance. In all cases, 
there is an initial transient period, 15 to 20 cycles long, where the TBF increases by about a third before 
settling into an indefinite steady state with constant frequency. We subsequently re-inspected behavioral 
sequences, and fou nd varying degrees of wind-up in some, but not all, experimentally recorded swim bouts 
l)Borla et all 1200^ . It is unknown whether the biological and simulational wind-up behaviors are related, 
but the phenomenon lasts much longer in the model than in the lab. 

To evaluate the range of oscillator or tail-beat frequencies that might be generated with this simple 
model, we tested different combinations of synaptic strengths (conductances). Varying the strengths of 
the NMDA and glycinergic synapses was found to alter the steady-state oscillator frequency (Fig. 3): 
increasing the NMDA conductance increased TBF, whereas increasing the glycinergic conductance de- 
creased it. The AMPA conductance can also affect TBF, but to a lesser extent. Different combinations 
of conductances could give rise to the entire range of TBFs observed in different larval swim patterns, 
from the slow swim (25 to 40 Hz) to burst swim (45 to 75 Hz: .Budick and O'Malkv. 2000) . This is just 
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Fig. 3. A density plot showing the steady-state frequency of a single oscillator as we modify the synaptic conductances 
9NMDA a-iid 9lNHi fixing (;ampa = 10~*. All conductances are in /iS, and all frequencies are in Hz. The missing regions in 
the upper-right and lower-left corners represent non-oscillatory states. 
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Fig. 4. Rastergram for a chain of 25 segmental oscillators. Each horizontal line corresponds to a single oscillator; ticks drawn 
below each line are for action potentials on the left, and above for those on the right. Parameters: (?ampa = 10~* /iS, 

9NMDA = 6 X 10-" /iS, giNH = 10-2 ^S, 3doscending = 10-^ fjS. 

one potential means of varying oscillator frequency (or TBF), implemented in a reduced system (i.e. a 
two-cell model), but it illustrates how a minimal model of larval zebrafish spinal cord, with just a few 
essential conductances, can give rise to a range of outputs relevant to the larval behaviors. In this model, 
each crossed inhibitory signal results in one post-inhibitory rebound firing of a single action potential, but 
to more completely capture slow and burst swims we will need to incorporate mechanisms that regulate 
the strength of output of the motoneuron pools. 

We next extended the model to a 25-segment chain of oscillators, creating a kind of artificial spinal 
cord that should facilitate quantitative analyses of the influences of descending signals on spinal network 
activity. The firing pattern of the model spinal cord is illustrated in Fig. 4. The firing of segment num- 
ber 1 shows an alternating left-right pattern over the duration of the simulation. In successively more 
caudal segments, the firing is delayed, matching the general pattern of undulatory swimming. As in the 
single segment model, each segment along the chain showed a frequency wind-up period. In this model a 
fixed intersegmental time-delay (nominally a synaptic delay) was used to establish the phase relationship 
between segments, but a more realistic mod el might include ascen ding an d descendi ng connections of 
varying strengths and lengths, as in lamprey ijKotaleski et alJll999l and see 

The kinematic model of the larval trunk is a complementary tool for exploring theories of descending 
motor control. Preliminary work (;Hill et al.. . .200,11 assumed that the fish was equally flexible from trunk 
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Fig. 5. Comparison of a 7-day old larval zobrafish with our kinematic model during a slow swim (left) and a J-turn (right). 

to tail {W{x) = 1 in Eq. 2), with partial success. By accounting for the increased flexibility of the tail 
(e.g. by making W{x) linear), however, we can better simulate the trunk kinematics observed during 
slow swims and J-turns (Fig. 5). J-turns are uniqu e locomotive maneuvers that contribute to the larval 
prey-capture behavior ijBorla and O'Malievl l2002(l . They require repetitive, asymmetric and far-caudal 
contractions of axial musculature. This might in principle be achieved by sending an asymmetric excitatory 
brainstem signal exclusively to far-caudal spinal cord; however, a survey of the spi nal outputs of zebrafish 
reticu lospinal neurons revealed no neurons with the requisite arborization pattern ijGahtan and O'Malievl 
|2003|) . Alternatively, neurons that selectively arborize in rostral spinal cord (which were observed) might 
be activated bilaterally to stiffen the rostral musculature. In conjunction with such signals, other neurons 
that arborize along the entire length of spinal cord could then generate far-caudal contraction. In our 
model, this is implemented as a bilateral "inhibitory" rostral signal, but this is kinematically equivalent 
to stiffening the rostral end of the larva by bilateral excitation of rostral motoneurons. 

Future Directions 

As we extend the neural model's capabilities to produce dynamically varying bend amplitudes, we 
hope to produce an increasingly realistic "artificial spinal cord" that can be used to test ideas of how the 
larval locomotive repertoire is generated. The kinematic model is useful here because it provides a first 
approximation of how complex neural model outputs might affect axial kinematics. Further extensions 
of this model should incorporate known constraints of the larval CNS, such as the diversity of brainstem 
and spinal systems involved in swimming and turning behaviors (jO'Mallev et al.lll99 6: Hale ct al., 2001; 
Ritter et al.. 2001). and the wide distribution of activity during escape behaviors (Bosch ct_al,,, 2001 ; 

ahtan et al. . l20()2) . By incorporating such constraints, the combined neurokinematic model should 
become increasingly useful in generating experimentally-testable hypotheses of descending and spinal 
control of vertebrate locomotion. 
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